The Nuclear Matter Symmetry Energy at 0.03 < p/po < 0.2 
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Measurements of the density dependence of the Free symmetry energy in low density clustered 
matter have been extended using the NIMROD multi-detector at Texas A&M University. Thermal 
coalescence models were employed to extract densities, p, and temperatures, T, for evolving systems 
formed in collisions of 47 A MeV ''°Ar + ^^^Sn and ^^Zn + "^Sn , ^^''Sn. Densities of 

0.03 < p/po < 0.2 and temperatures in the range 5 to 10 MeV have been sampled. The Free 
symmetry energy coefficients are found to be in good agreement with values calculated using a 
quantum statistical model. Values of the corresponding symmetry energy coefficient are derived 
from the data using entropies derived from the model. 

PACS numbers: 25.70.Pq 

Keywords: Intermediate Heavy ion reactions, chemical equilibrium, neutron and proton chemical potential, 
quantum statisitcal model calculations 



I. INTRODUCTION 

The symmetry energy of nuclear matter is a fundamen- 
tal ingredient in the investigation of nuclear and astro- 
physical phenomena. The symmetry energy character- 
izes the dependence of the nuclear binding energy on the 
asymmetry i5 = {N — Z)/A where Z and N are the pro- 
ton and neutron numbers, and A = N + Z . As a general 
representation of the symmetry energy coefficient we use 
the definition 

E,y^{p, T) = i {E{p, 1, T) + E{p, -1,T))~ E{p, 0, T) 

(1) 

where E{p, 6, T) is the energy per nuclcon of nuclear 
matter with density p, asymmetry 5, and temperature 
T . If a quadratic dependence is assumed this definition 
becomes identical to the frequently used alternative of 
taking the second derivative of E{p,6,T) with respect 
to the asymmetry 6 = 0. Our empirical knowledge of 
the symmetry energy near the saturation density, po, is 
based primarily on the binding energies of nuclei. The 
Bethe-Weizsacker mass formula leads to values of about 
^^sym(po: 0) = 28 — 34 McV for the symmetry energy at 
zero temperature and saturation density po = 0.16 fm"-^, 
if surface asymmetry effects are properly taken into ac- 



count [l|. In contrast to the value of £'sym(/5o,0), the 
variation of the symmetry energy with density and tem- 
perature is intensely debated. Many experimental and 
theoretical investigations have been performed to esti- 
mate the behavior of the symmetry energy as a function 
of p and T. Recent reviews are given by Li et al. Q and 
by Lattimcr and Lim 

In the Fermi energy domain symmetry energy effects 
have been investigated using judiciously chosen observ- 
ables from nuclear reactions [2|-|7|. In the theoretical in- 
vestigations quasiparticle approaches such as the Skyrme 
Hartree-Fock and relativistic mean field (RMF) models 
or Dirac-Brueckner Hartree-Fock (DBHF) calculations 
are used 0, 0: @- In such calculations the symmetry 
energy tends to zero in the low-density limit for uniform 
matter. However, in accordance with the mass action 
law, cluster formation dominates the structure of low- 
density symmetric matter at low temperatures. Thus at 
low density the symmetry energy changes mainly because 
additional binding is gained in symmetric matter due to 
formation of clusters and pasta structures [l3|. As a re- 
sult the symmetry energy in this low-temperature limit 
has to be equal to the binding energy per nucleon as- 
sociated with the strong interaction of the most bound 
nuclear cluster. A single-nucleon quasiparticle approach 
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cannot account for such structures. The correct low- 
density hmit can be recovered only if the formation of 
clusters is properly taken into account, as has previously 
been shown in [111 ] in the context of a virial expansion 
valid at very low densities and in Ref. [l2j. Since, at 
low density, the symmetry energy changes mainly be- 
cause additional binding is gained in the formation of 
clusters [HI, measurements of nucleon and light 

cluster emission from the participant matter which is pro- 
duced in near Fermi energy heavy ion collisions can be 
employed to probe the EOS at low density and moder- 
ate temperatures where clustering is important fl3. 14 1. 
Our previous data demonstrated a large degree of al- 
pha clustering for densities at and below 0.05 times 
normal nuclear density, (0.16 nucleons/fm"^) and tem- 
peratures of 4 to 10 MeV. Using these data we derived 
symmetry energy coefficients in low density nuclear mat- 



ter [13|, |lj| . The analysis employed the isoscaling tech- 
nique which compares yields for two systems with simi- 
lar temperatures but different N jZ ratios to determine 
the differences in chemical potentials and symmetry en- 
ergy [Hill. The NIMROD 47r muhi-dctcctor at Texas 
A & M University has now been used to extend our mea- 
surements to higher densities. Cluster production in col- 
lisions of 47 A McV 4°Ar with "2424gj^ ^^^^ M-^n with 
ii2,i24gj-^ studied. We report here determinations of 
symmetry energy coefficients at 0.03 < p/po < 0.2 and 
moderate temperatures. Our results for this expanded 
range of densities are in reasonable agreement with those 
of a quantum statistical model calculation incorporat- 
ing medium modifications of the cluster binding energies 



12, la 2i 



EXPERIMENTAL TECHNIQUES 

The experiments were performed using beams from 
the K500 Superconducting Cyclotron at the Texas A&M 
University Cyclotron Institute, incident on targets in the 
NIMROD detector. NIMROD consists of a 166 segment 
charged particle array set inside a neutron ball [2lj . The 
charged particle array is arranged in 12 rings of Si-Csl 
telescopes or single Csl detectors concentric around the 
beam axis. The Csl detectors are 1-10 cm thick Tl doped 
crystals read by photomultiplier tubes. A pulse shape 
discrimination method is employed to identify light parti- 
cles in the Csl detectors. For this experiment each of the 
forward rings included two segments having two Si de- 
tectors (150 and 500 /im thick) in front of the Csl detec- 
tors (super telescopes) and three having one Si detector 
(300 /im thick). Each super telescope was further divided 
into two sections. Neutron multiplicity was measured 
with the 47r neutron detector surrounding the charged 
particle array. This detector is a neutron calorimeter 
filled with gadolinium doped pscudocumene. Thermal- 
ization and capture of emitted neutrons in the ball leads 



to scintillation light which is observed with phototubes 
providing event by event determinations of neutron mul- 
tiplicity. Further details on the detection system, energy 
calibrations and neutron ball efficiency may be found in 



references [2^ l23| . The combined neutron and charged 



particle multiplicities were employed to select the most 
violent events for subsequent analysis. 



ANALYSIS 



The dynamics of the collision process allow us to probe 
the nature of the intermediate velocity "nucleon-nucleon" 

2? 



emission source 
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Measurement of emis- 
sion cross sections of nucleons and light clusters together 
with suitable application of a coalescence ansatz |24| pro- 
vides the means to probe the properties and evolution of 
the interaction region. The techniques used in the anal- 
ysis have been detailed in several previous publications 
13l Il4l . I22I |23| and are described briefly below. A no- 
table difference from Rcfs. 131 and \\M is the method 



of density extraction. This is discussed more extensively 
in the following. We emphasize that the event selection 
is on the more violent collisions. Cross section weighting 
favors mid-range impact parameters. 

An initial estimation of emission multiplicities at each 
stage of the reaction was made by fitting the observed 
light particle spectra assuming contributions from three 
sources, a projectile-like fragment (PLF) source, an inter- 
mediate velocity (IV) source, and a target-like fragment 
(TLF) source. A reasonable reproduction of the observed 
spectra is achieved. Except for the most forward detec- 
tor rings the data are dominated by particles associated 
with the IV and TLF sources. The IV source velocities 
are very close to 50% of the beam velocity as seen in many 
other studies 
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The observed spectral slopes re- 
flect the evolution dynamics of the source, not its internal 
temperature [2^, [2^. For further analysis, this IV source 
is most easily sampled at the intermediate angles where 
contributions from the other sources are minimized. For 
the analysis of the evolution of the source we have se- 
lected the data in ring 9 of the NIMROD detector. This 
ring covered an angular range in the laboratory of 38° to 
52°. An inspection of invariant velocity plots constructed 
for each ejectile and each system, as well as of the results 
of the three-source flt analyses, indicates that this selec- 
tion of angular range minimizes contributions from sec- 
ondary evaporative decay of projectile like or target like 
sources. We treat the IV source as a nascent fireball cre- 
ated in the participant interaction zone. The expansion 
and cooling of this zone leads to a correlated evolution of 
density and temperature which we probe using particle 
and cluster observables, yield, energy and angle. 
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Temperature 



As in some previous work 13, 14] we have employed 
double isotope yield ratios (27l 28 1 to characterize the 
temperature at a particular emission time. For parti- 
cles emitted from a single source of temperature, T, and 
having a volume Maxwellian spectrum the HHe double 
isotope yield ratio evaluated for particles of equal Usmf is 
•y/9/8 times the ratio derived from either the integrated 
particle yields or the yields at a given energy above the 
barrier. 



14.3 



HHg 



ln(V9/8(1.59 



(2) 



If Y represents a cluster yield, Rv^^^t = 
y(2H)y(4He)/r(3H)y(3He) for clusters with the 
same surface velocity. The constants 14.3 and 1.59 
reflect binding energy, spin, masses and mass differences 
of the ejectiles. Eq. ([2|) differs from the usual formulation 
only by the factor of y^9/8 appearing in the logarithm 
term in the denominator. 

Model studies [2^ comparing Albergo model 27 1 tem- 
peratures and densities to the known input values have 
shown the double isotope ratio temperatures to be rela- 
tively robust in this density range. However the densities 
extracted using the Albergo relationships are useful only 
at the very lowest densities. Consequently, in this study 
we have employed a different means of densit y ex traction 
- the thermal coalescence model of Mekjian 
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Density 

To determine the coalescence parameter Pq , the radius 
in momentum space, from our data we have followed the 
Coulomb corrected coalescence model formalism of Awes 
et al. (soj and previously employed by us in Ref. psf . 
In the laboratory frame the derived relationship between 
the observed cluster and proton differential cross sections 
is 
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dEAdn 



R 



N 1 

"PNIZIA 
d'^N{l,0,E) 



47rP3 



3[2m3(£; 

A 



dEdn 



A-l 



(3) 



where the double differential multiplicity for a cluster of 
mass number A containing Z protons and N neutrons 
and having a Coulomb-corrected energy Ea, is related 
to the proton double differential multiplicity at the same 
Coulomb corrected energy per nucleon, E — Ec, where 
Ec is the Coulomb barrier for proton emission. i?„p is 
the neutron to proton ratio. 
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FIG. 1: (Color online) Coalescence parameters, Po as a func- 
tion of surface velocity in the intermediate velocity source 
frame. Reaction: 47 A MeV ''°Ar + "^Sn 



Application of the coalescence model requires knowl- 
edge of cluster, neutron and proton differential cross sec- 
tions with proper absolute normalizations. In this work 
absolute measured multiplicities for the selected violent 
events are employed. The neutron spectra are not mea- 
sured. However, since within the framework of the coa- 
lescence model the yield ratios of two isotopes that differ 
by one neutron are determined by their binding energies 
and the n/p ratio in the coalescence volume, we have 
used the observed triton to '^He yield ratio to derive the 
n/p ratio used in this analysis. Because our goal was to 
derive information on the density and temperature evolu- 
tion of the emitting system, our analysis was not limited 
to determining an average Pp value. Instead, as in our 
previous studies [H-liillii,!!!, results for d, t, ^He, and 
^Hc were derived as a function of the surface velocity in 
the intermediate velocity source frame, i.e., the velocity 
of the emerging particle at the nuclear surface, prior to 
Coulomb acceleration. Results for the reaction of ^°Ar 
+ ^^^Sn are presented in Fig. [TJ Values of Pq for the 
other reactions are very similar. The low velocity cut off 
is imposed to avoid significant contamination from par- 
ticles evaporated from the TLF source. The velocity of 
6.6 cm/ns in the IV frame corresponds to the beam ve- 
locity and the rise in Pq seen in the Figure indicates that 
these high velocity particles may not be amenable to an 
equilibrium treatment. 

In the Mekjian model thermal and chemical equilib- 
rium determines coalescence yields of all species. Under 
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FIG. 2: (Color online) Coalescence model volumes as a func- 
tion of velocity in the intermediate velocity source frame. Re- 



action: 47 A MeV ''°Ar + "^Sn 



these assumptions there is a direct relationship between 
the derived radius in momentum space and the volume 
of the emitting system. In terms of the Pq derived from 
Eq. ([3]) and assuming a spherical source, 



V = 



4^ 



ZINIA^ 
2^ 



(2s + l)eT^ 



(4) 



where h is Plancks constant, and Z, N, and A are the 
same as in Eq. (|3]), i?o is the ground state binding energy, 
s the spin of the emitted cluster, and T is the tempera- 
ture. Thus the volume can be derived from the observed 
Pq and temperature values assuming a spherical source. 
We note that this volume is a free volume. From the rele- 
vant Pq values we then determined volumes using Eq. (j4|) . 
See Fig. [5] where, again, results for the reaction Ar + 
^^^Sn are presented. 

A comparison of these volumes indicates good agree- 
ment for t, ■^Hc and "'He. Volumes derived from the 
deuteron data are typically somewhat larger. This larger 
apparent volume for deuterons appears to reflect the 
fragility of the deuteron and its survival probability once 
formed [sH. For this reason we have used average vol- 
umes derived from the A > 3 clusters to calculate the 
densities. Given that mass is removed from the system 
during the evolution, we determined the relevant masses 
for each volume by assuming that the initial mass of the 
source was that obtained from the source fitting analy- 
sis and then deriving the mass remaining at a given Vsmf 
from the observed energy spectra. This is also an averag- 



ing process and ignores fluctuations. Once these masses 
were known they were used to determine an excluded vol- 
ume for the particles. Addition of this excluded volume 
to the free volume produced the total volumes needed 
for the density calculations. These were determined by 
dividing the remaining masses by the total volume. This 
was done as a function of Wsurf- This excluded volume 
correction is small for low Vsmf but increases with in- 
creasing fsurf- Densities are those of total number of 
nucleons, including the nucleons bound in clusters, per 
fm'^. 



RESULTS 



Temperatures and Densities 



Inspection of the results for the four different systems 
studied revealed that, as a function of surface velocity, 
the temperatures, densities and equilibrium constants for 
all systems are the same within statistical uncertainties. 
This indicates that the time evolution of the systems 
studied is very similar (23l . Is^ . Therefore we have com- 
bined values from all systems to determine the tempera- 
tures and densities reported in this paper. 

We present in Fig.[3]the experimentally derived density 
and temperature evolution for the IV source. Estimated 
errors on the temperatures are 10% at low density evolv- 
ing to 15% at the higher densities. In deriving our results 
we have invoked the correlation between average surface 
velocity and emission time for the early emitted parti- 
cles on which we focus our attention. This is seen clearly 
in transport model calculations [s^. It is certainly true 
that density and temperatures corresponding to the dif- 
ferent surface velocity bins are weighted averages over 
the underlying distributions. There will be fluctuations 
in T and p present at the time of particle emission. We 
assume these to be symmetric about the most probable 
value. There will then be additional fluctuations induced 
by mixing. Using the AMD model of Ono, and assum- 
ing at each time a single (time-decreasing) temperature 
we estimate that, for the higher surface velocities, the 
weighted average temperatures are < 10% lower than the 
input temperatures. This effect decreases at lower sur- 
face velocities where the rate of decrease of temperature 
with velocity is less. These estimates are included in our 
estimates of the temperature uncertainty. The error in 
the derivation of the density arises from the uncertainty 
on the volume which is dominated by the uncertainty in 
temperature and the uncertainty in source mass derived 
from source fitting to complex spectra. The estimated 
errors on the densities arc ±17%. 
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FIG. 3: (Color online) Temperatures and densities sampled 
by cluster emission from the expanding IV source. 



Isoscaling 

Considering only bulk properties of homogeneous nu- 
clear matter, and the standard thermodynamic relations 
that relate the Free energy, the Internal energy, and the 
entropy, as a function of p, 5, T we can expand with re- 
spect to 5 to define the corresponding expressions for 
F,y^{p,T) and 5,y„,(p,T) = -dF,y^{p,T)/dT. This is 
analogous to the expressions given above for Esym{p,T) 

The isoscaling technique has been used to derive in- 
formation on Fsym , the symmetry Free energy of the sys- 
tem. The isoscaling technique consists of measuring clus- 
ter yield ratios for two different excited systems (denoted 
by the index 1 and 2), having the same temperature (and 
density) and the same atomic number but differing in the 
Z/A ratio. For such systems liquid drop model binding 
energy terms other than those related to the symmetry 
energy coefficient essentially cancel and one can write 



^2 _ (^g((M2(n)-Ml("))A^+(A'2(p)-Ml(p))^)/T 



Yi 



Ce 



aN+pZ 



(5) 



where C is a constant and p{n) and p{p) are the neutron 
and proton chemical potentials. The isoscaling parame- 
ters a = {p.2{n)-pi{n))/T and /? = {p^ip) - Pi{p)) / T , 
representing the difference in chemical potential between 
the two systems, may be extracted from suitable plots of 
yield ratios. Either parameter may then be related to the 
symmetry Free energy, i^sym- We take the a parameter. 



which is expected to be less sensitive to residual Coulomb 
effects. Addressing specifically the symmetry Free energy 
and adopting the usual convention that system 2 is richer 
in neutrons than system 1, one can write 



4K, 



A2 



IT. 



(6) 



where Z is the atomic number and A is the mass number 
of the emitter [H, 17, 18|. Thus, i^sym may be derived di- 
rectly from determinations of system temperatures, Z j A 
ratios, and isoscaling parameters. We emphasize that the 
present analysis is carried out for light species character- 
istic of the nuclear gas rather than, as in most previous 
analyses, for the intermediate mass fragments thought to 
be characteristic of the nuclear liquid. In this work we 
employ Eq. ^ with experimentally determined isoscal- 
ing parameters, a, temperatures, T, and ZfA ratios to 
determine the symmetry Free energy coefficient, i^gym- 

Fig. U (a-e) presents the results of our isoscaling analy- 
sis. The experimentally derived evolution of density, tem- 
perature, \{Z\I Aif' ~{Z2l A2)^\^ the isoscaling parameter 
a, and the Free symmetry energy coefficient are presented 
as a function of surface velocity in the IV frame. The 
correlations between these parameters become apparent 
in this figure. For the determination of the Free sym- 
metry energies we have restricted ourselves to isoscaling 
results leading to alpha parameters with less than 15 % 
uncertainties. This restricts the temperature and density 
range for which the Free symmetry energies are reported 
in figure 4. 

We see that T and p increase with increasing fgurf ■ At 
the same time there is a decline in the isoscaling coeffi- 
cient a. Combined, the terms of Eq. ([5]) lead to values of 
^sym that rise gently with 

'^surf ■ 



Free Symmetry Energy 

In Figure 5 the Free symmetry energy coefficients de- 
termined from the isoscaling analysis are plotted. Also 
shown are calculated Free symmetry energy coefficients 
obtained from a quantum statistical approach that in- 
cludes medium modifications of the cluster binding ener- 
gies 12, 15, 3, 3, i^l- A few-body Schrodinger equation 



has been solved that contains single-particle self-energy 
shifts as well as Pauli blocking factors. For given T and 
Pp,Pm the nucleon density is calculated. To obtain a 
thermodynamic potential, the Free energy as function of 
T, n„,7T.p is evaluated by integration. The further ther- 
modynamic quantities are determined from the Free en- 
ergy in a consistent way. 

Isotherms of the Free symmetry energy at T = 4, 6, 8 
and 10 MeV are also shown in Fig. 5. Note that the ex- 
perimental data correspond to changing temperatures as 
seen in Fig. 4. Within the error bars, the general behav- 
ior of the experimental data is fairly well reproduced by 
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FIG. 4: Parameters of the analysis, (a) Symmetry Free energy 
coefficient, (b) Temperature, (c) Isoscaling parameter, a, (d) 
(f^)' - (If)', (e) Density 



FIG. 6: (Color online) Symmetry entropy vs density. Symbols 
show experimental points and lines show calculations for T=4, 
6, 8 and fO MeV. 
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FIG. 5: (Color online) Free symmetry energy vs. density. 
Symbols show experimental points. Lines are calculated val- 
ues for T =4, 6, 8 and 10 MeV fTl, [TE, [3, [ll, [li 



the calculations. The data are somewhat below the cal- 
culated values, in particular in the intermediate region. 

To derive the values of the Internal symmetry energy 
from the Free symmetry energy, the symmetry entropy 
must be known. The calculations of the symmetry en- 
tropies have been performed using the model of Typel et 
al. [12[ where the Free energy, the Internal energy and the 
entropy have been given for symmetric matter. For pure 
neutron matter, the RMF approach can be used because 
there is no cluster formation. 

We calculated the entropy within the QS approach by 
differentiating the Free energy with respect to the tem- 
perature. In Figure 6 we present the symmetry entropy, 
i.e. the difference between the entroy of neutron matter 
and that of symmetric matter, for different temperatures 
using the recent quantum statistical calculations of mo- 
mentum dependent energy shifts of light clusters [16 1. 

In contrast to the mixing entropy that leads to a larger 
entropy for uncorrelated symmetric matter in comparison 
with pure neutron matter, the formation of correlations, 
in particular clusters, will reduce the entropy in symmet- 



ric matter, see also Fig. 9 of Ref. [13]. For parameter 
values for which the yields of free nucleons in symmet- 
ric matter are small, the symmetry entropy may become 
positive. The fraction of nucleons bound in clusters can 
decrease due to increasing temperature or the dissolu- 
tion of bound states at high densities due to the Pauli 
blocking. Then, the symmetric matter recovers its larger 
entropy so that the symmetry entropy becomes negative. 
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FIG. 7: (Color online) Symmetry energy vs density. Symbols 
show experimental points and lines show calculations for T=4, 
6, 8 and 10 MeV. 



The corresponding results for the symmetry energy are 
shown in Fig. 7. Experimental as well as calculated val- 
ues indicate a large value of the symmetry energy also 
in the low-density region since the alpha-fraction in sym- 
metric matter becomes large at low temperatures, and 
the value of the binding energy per nucleon determines 
the Internal energy. This effect in not included in quasi- 
particle approaches such as standard relativistic mean 
field calculations because few-body correlations are ne- 
glected. 



SUMMARY AND CONCLUSIONS 

The NIMROD multi-detector at Texas A&M Univer- 
sity has been employed to extend our earlier measure- 
ments of symmetry energy coefficients in low density clus- 
tered nuclear matter. Yields of light particles produced 
in the collisions of 47 A MeV ^^Ar with ^^^Sn ,i2'^Sn 
and ^^Zn with ^^^Sn, ^^''Sn were employed in Thermal 
coalescence model analyses to derive densities and tem- 
peratures of the evolving emitting systems. Isoscaling 
analyses were used to determine the Free symmetry en- 
ergies of these systems. Comparisons of the experimental 
values are made with those of calculations made using a 
model which incorporates medium modifications of clus- 



ter binding energies. The model calculated symmetry 
entropies have been used together with the experimen- 
tal Free symmetry energies to derive symmetry energies 
of nuclear matter at densities of 0.03 < p/po ^ 0.2 and 
temperatures in the range 5 to 10 MeV. 
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